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The measurement of diffractively scattered protons in the ATLAS Forward Physics 
detector system placed 220 m away from the ATLAS interaction point is studied. A 
parameterisation of the scattered proton transport through the LHC magnet lattice is 
presented. The proton energy unfolding and its impact on the centrally produced scalar 
particle mass resolution are discussed. 

1. Introduction 

Diffractive dissociation is one of the processes that can be studied at the LHC. Diffrac- 
tive physics is strong interaction physics involving no exchange of quantum numbers other 
than those of the vacuum. In experiment this leads to an obvious triggering scheme relying 
on the rapidity gap method. However, one has to keep in mind that the gap has a certain 
survival probability that depends on the interaction type and the centre of mass energy. 
In addition, also the diffractively scattered protons can be tagged. Usually, the protons 
scatter at small angles and in the collider environment they stay in the beam pipe and 
travel through the magnet lattice of the machine. A possible way to measure parameters 
of the diffractively scattered proton trajectory is to use detectors placed inside the beam 
pipe, for example by means of a roman pot station or a movable beam pipe technique. 
The ATLAS Collaboration plans to have proton tagging stations placed symmetrically 
with respect to the Interaction Point (IP) at the distances of 220 m and 420 m (AFP220 
and AFP420) and 240 m (ALFA). The ALFA (Absolute Luminosity For Atlas) [Q] sta- 
tions at 240 m will be devoted to the absolute luminosity measurement of the LHC at 
the ATLAS IP. This measurement will rely on the detection of the elastically scattered 
protons. The AFP (Atlas Forward Physics) [ |2] stations will be used for diffractive and 
77 physics. 
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2. Experimental Environment 

Below, only the aspects concerning the AFP220 detectors are discussed. In the mea- 
surement, the machine magnets play the role of magnetic spectrometer. Therefore, two 
detector stations are to be placed around 220 m, namely one at 216 and one at 224 m. 
Each station will be equipped with position sensitive and triggering detectors, horizon- 
tally inserted into the LHC beam pipe. The position sensitive detectors will consist of 
10 layers of the silicon 3D detectors [[3] to measure the scattered proton trajectory. The 
position measurement resolutions are assumed to be <j x = 10 /xm and <y y = 40 /im in the 
horizontal and vertical direction, respectively. This will allow the measurement of the 
particle trajectory position and direction. Additionally, the stations will contain very fast 
timing detectors with picosecond resolution. They will measure the scattered proton time 
of flight and indicate the interaction vertex longitudinal coordinate, z, with resolution 
of the order of several millimeters. In the following calculations, a reference frame with 
the x-axis pointing towards the accelerator centre, the y-axis pointing upwards and the 
z-axis along one of the beams was used. 

The scattered proton can be described at the interaction point in several equivalent 
ways, each useful in a different case: (p x ,Py,Pz), (E,x' ,y' ), or (E,p T ), 

E = A/m 2 + p 2 x = — y = — 

Pz Pz 

where E and m are the proton energy and mass, p = {p x ,Py,Pz) is the proton momentum 
and px = {Px,Py) denotes the proton transverse momentum. Useful variables are the 
proton energy loss AE = E — E (E is the incident beam energy), the reduced proton 
energy loss, £ = AE/Eq, and the four- momentum transfer between the incident and 
scattered proton, t. 

There are several programs on the market calculating proton trajectories through the 
magnets. In the following the FPTrack program [ H| was used. This program computes 
the positions of particles using the LHC optics files which describe the magnetic fields, 
positions and apertures of the LHC lattice. These files were produced with help of a 
principal beam transport program Mad-X [[5] by the LHC optics group [[6]. 

Table 1 

The LHC beam and the crossing region parameters at the ATLAS IP. 



Parameter 


Beam 


Crossing region 


°x 


16.6 iim. 


11.7 iim. 


a yo 


16.6 ixm. 


11.7 /im 


<**o 


75 mm 


53 mm 


°x> 


30.2 /xrad 




a y'o 


30.2 /zrad 




°E 


0.77 GeV 
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The interaction vertex is described by its coordinates (xq, yo, Zq). These coordinates 
have Gaussian distributions with zero average values and dispersions: a XQ , a yo , a zo , re- 
spectively. In the simulation, the beam particle energy and its momentum direction 
were generated according to Gaussian distributions with appropriate means and the dis- 
persions: a Eo for the energy and a x ' Q ,a y ' o for the momentum in (x,z) and (y,z) planes, 
correspondingly. Values of these parameters for the nominal 7 TeV beam energy and 
standard LHC optics are listed in Table [TJ The calculations were performed for both 
beams: beaml that performs the clockwise motion and beaml which does the counter 

cJnnkwisp rotation 



a) 



beam 1 



b) 




beam 2 



-2-1012 -2-1012 

x [mm] x [mm] 

Figure 1. The LHC beam profiles at 216 meters from ATLAS IP for beaml (a) and beam2 
(b), for the 7 TeV LHC optics. 



Figure [T] shows the LHC beam profiles in the (x,y) plane at 216 m away from the 
ATLAS IP obtained with the FTPtrack program for both beams and the 7 TeV LHC 
optics. As can be observed, the beams are much wider in the vertical direction than in 
the horizontal one. 

Table [2] summarizes the beams spreads following from a two dimensional Gaussian fit to 
the beam profiles shown in Fig. [TJ The (10 - 15)a beam envelope gives a natural limit for 
the distance between the detector frame and the beam centre. In the horizontal direction 
this corresponds to about 1-2 millimeters. Obviously, this distance plays a crucial role 
for the diffractively scattered proton detection and hence, for the experimental apparatus 
acceptance. 

It is important to see in which range of the energy, E, and the transverse momentum, 
Pt = \Pt\, the detector can measure protons. The geometric acceptance for fixed E and 
Pt values was defined as the ratio of the number of protons that crossed the detector to 
the total number of scattered protons with a given E and pt- Only the effects of the beam 
pipe aperture and the distance between the detector and the beam centre were taken into 
account. 

Figure [2] depicts the geometric acceptance as a function of E and pr for beaml for the 
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Table 2 

The LHC beam spreads from a Gaussian fit to the beam profiles at 216 m. 



Parameter 
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Beam 2 
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Figure 2. The geometrical acceptance of the detector placed in the LHC beaml at 216 
m away from the IP as a function of the proton energy loss {AE) and its transverse 
momentum px for a 3 mm distance between the beam centre and the detector active 
edge. 



standard 7 TeV LHC optics. The acceptance is above 80% in the region limited by: 

200 < AE < 1000 [GeV], < p T < 2.5 [GeV/c] 
which corresponds to 

0.03 < f < 0.14, -6.5 < t < [GeV 2 /c 2 ]. 

When the geometrical acceptance requirement is lowered to 60%, this results in a wider 
range of the accepted proton energies and transverse momenta. The range enlargement 
is seen (c.f. Fig. ED for 200 < AE < 600 GeV and p T < 3 GeV/c. This gives the limits: 

0.03 < £ < 0.14, -10 < t < [GeV 2 /c 2 ]. 

The presence of two detector stations on each side of the ATLAS detector allows the 
measurement of the proton trajectory elevation angles, x ' and y ', in the (x,z) and (y,z) 



Transport Simulation and Diffractive Event Reconstruction at the LHC 



5 



planes, respectively. From the FPTrack calculations, it follows that the position and slope 
of the trajectory at the detector in one transverse direction is independent of those in the 
other direction, i.e. x and x' values do not depend on y and y' Q and vice versa. This is 
a reflection of a negligible role of the sextupole and higher order magnetic fields in the 
standard LHC optics between the ATLAS IP and the AFP220 stations. 
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Figure 3. The x-direction chromaticity plots for the LHC beaml (a) and beam2 (b) for the 
standard 7 TeV LHC optics. The lines of constant energy correspond to 7000, 6825, 6650, 
6475, 6300 GeV from left to right, respectively. The angles were changes from -400 /zrad 
to 400 yurad (from top to bottom). 



To illustrate how the proton trajectory positions and slopes measured by the detectors 
depend on the proton energy and its trajectory slopes at the IP the chromaticity plots 
were prepared. The plots shown in Figure [3] were devised by plotting in the (x,x r ) plane 
the lines corresponding to the constant E and x' at the IP. 

The chromaticity plots indicate few things. Firstly, there is a non-negligible difference 
between both beams. Therefore, properties of both have to be studied. Secondly, the grids 
created by the energy and angle iso-lines do not fold. Hence, it is possible to obtain the 
energy and the transverse momentum of a proton from the measurements of the proton 
trajectory in both stations. In particular, assuming a fixed interaction vertex position (no 
smearing) the energy of a proton at the IP can be deducted solely from the measured x 
and x' values. 

3. Transport Parameterisation 

In order to unfold the proton energy from the detector measurements a parameterisation 
of the FPTrack transport calculations was prepared. The aim was to describe the FPTrack 
results analytically. It was requested that: 

• the parameterisation has a simple functional form, 

• the parameterisation precision has an accuracy which is better than the assumed 
detector spatial resolutions. 
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To find the parameterisation form, events uniformly distributed over the (E,x' Q ,y' Q , 
x ,y ,z ) space were generated. Subsequently, these events were used in the FPTrack 
transport calculations. The transport results were the input data to the parameterisation 
search procedure. It was found that the following parameterisation fulfill the requirements 
outlined above well: 

a = A a + a' B a + a C a + a' z D a + z F ai (1) 

a ' = A sa + a' B sa + a C sa + a' z D sa + z F sa , (2) 

where a = {x,y}, sa denotes the slope either in x or y direction and all the capitalised 
symbols are polynomials of energy, i.e.: 



The values of all the coefficients were found by fitting the formulae to the FTPtrack 
calculations for simulated events. 

The accuracy of the method was estimated by plotting the difference between the 
value given by the parameterisation and that given by the FPTrack calculation. The 
accuracy of the position parameterisation was found to be of the order of a microme- 
ter which is 10 times less than the assumed detector resolution in the horizontal plane. 
The difference between FPTrack and the parameterisations for the trajectory angles was 
found to be limited to about 50 nanoradians. One has to remember that the average 
multiple Coulomb scattering angle was estimated to be about 500 nrad. An example of 
the parameterisation accuracy is presented in Fig. |H In this figure the distributions of 
Ay 

— yparam y FPTrack 

and Ay' = y' param — y FPTrack are shown for single diffractive events 
generated with PYTHIA [[8]. The accuracy estimations for these quantities are displayed 
since they represent the worst precision cases. Nevertheless, the results are well confined 
within the ranges given by the detector resolutions. This confirms the parameterisation 
quality. 

One should note that the procedure outlined above can be easily repeated. In par- 
ticular, it can be applied to the files describing the actual LHC collision optics used for 
experimental runs. 

4. Event Reconstruction 

Since there is a correlation between the proton momentum and the measured position of 
the proton at the AFP220, the reconstruction of the proton properties at the interaction 
vertex from the measured coordinates of the proton trajectory at the AFP220 is possible. 

A proton at the interaction vertex is described by six independent variables: E, x' , 
y' , Xo, yo and z . The detectors deliver two pairs of transverse coordinates separated in 
longitudinal direction by a fixed distance of about 8 meters. In general, the unfolding 
problem is an ill-stated one. In the present case it requests the inversion of the 6 to 
4 mapping and its solution is only possible with help of additional assumptions. The 



A a = ai 0) + a^E + a^E 2 + a^E 3 + 



(3) 



Q = c (0) , (1) E _i_ c (2) E 2 , (3) E 3 



(4) 




Ay [/mi] Ay' [nrad] 

Figure 4. The parameterisation accuracy estimation examples (see text). Pictures show 
the uncertainty on the y (a) and y' (b) parameterisations of the beaml transport. 

simplest one is the assumption of a fixed position of the interaction vertex. In the following 
the positions x = y = and z = or z = 216 — cr were chosen, where r is the proton 
time of flight. 

To make the measurement simulation as close to reality as possible the detector ef- 
fects were taken into account. Protons traversing the detector station undergo multiple 
Coulomb scattering in the frame and the detector materials. The simulation of the proton 
trajectory position measurement also takes into account the assumed detector resolutions. 
These effects may lead to a considerable change of the proton trajectory parameters (slopes 
and positions). Hadron interactions in the detector or its frame were neglected as they 
are not important for the present study. 

A simple and fast method of the proton energy unfolding from the detector measure- 
ment is proposed. This method uses the assumption that the values actually measured 
are equal to those delivered by the parameterisation. This allows to calculate x' from 
eqs. (jTJ and (TSJ). Since both equations are considered for the same particle, they should 
give equal values. Hence, after simple algebra one gets: 

(x - A x - F x z - x C x ) ■ (B sx + z D sx ) = 

= (x' - A sx - F sx z - x C sx ) ■ (B x + z D x ) (5) 

where all capitalised symbols are described by eqs. [3] and EE 

The solution of the above equation is equivalent to finding the zero of the function f(E) 
given below: 

f(E) = (x- A x - F x z - x C x ) ■ (B sx + z D sx )- 

(x' - A sx - F sx z - x C sx ) ■ (B x + z Q D x ). (6) 

It was observed that for obtained parameterisation and simulated events the function 
f(E) has only one zero. Therefore, the equation 



f(E) = 
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can be easily solved numerically using for example the bisection method [ [7] . 

The energy unfolding procedure was tested using the same PYTHIA generated data 
sample. The proton energy was reconstructed with the help of the different additional 
assumptions listed below: 

• the "measured" trajectory coordinates were smeared according to the detector res- 
olution, 

• the interaction vertex transverse position was exactly known, 

• the interaction vertex longitudinal position was exactly known. 

The results are presented in Figure [5j The energy reconstruction resolution (the thick 
solid line) decreases from 9 GeV for 6000 GeV protons to about 3 GeV for 7000 GeV pro- 
tons. It is dominated by the detector spatial resolution which influence, marked with the 
thick dashed line, decreases with proton energy from about 7 GeV to about 1 GeV within 
considered energy range. Also, the impact of the multiple Coulomb scattering (dotted 
line) gets smaller with increasing proton energy. Its contribution to the resolution is about 
2.5 GeV at the maximum. For proton energies greater than 6800 GeV the uncertainty 
on the interaction vertex position in the transverse plane (the thick dash-dotted line) 



all effects B field variation 

detector resolution multiple scattering 

— ■ — xO unknown zO unknown 




6000 6200 6400 6600 6800 7000 

Proton energy [GeV] 

Figure 5. The proton energy reconstruction resolution for beaml as a function of its 
energy. The overall resolution is marked with the thick solid line, the influence of: the 
detector spatial resolution - the thick dashed line, the vertex position in the transverse 
plane - the thick dash-dotted line, the multiple Coulomb scattering - the dotted line, the 
vertex position along the beam axis - the dash-dotted line, the magnetic field variation - 
the solid line. 



Transport Simulation and Diffractive Event Reconstruction at the LHC 



9 



dominates the energy reconstruction resolution. The influence of the interaction vertex 
position along the longitudinal axis (the dash-dotted line) is small in the whole energy 
range discussed. 

Since the scattered protons can traverse the whole beam pipe volume the influence 
of the possible imperfections of the magnetic fields was studied. The magnetic fields of 
the lattice were varied by ±l%o of their nominal values. It should be pointed out that 
assumed variation is about a factor of 10 larger that the machine accepted and about 50 
times larger than the measured values [[9] of the higher multipoles at the reference radius 
of 17 mm away from the beam pipe centre. Variation of the magnetic field values gives a 
small contribution to the energy reconstruction resolution and for 6000 GeV protons it is 
about 0.3 GeV and increases to 0.7 GeV at 7000 GeV. This contribution is marked with 
the solid line in Fig. [51 The other effect of the variation of magnetic fields is the offset of 
the scattered proton reconstructed energy. This offset, on the absolute value, decreases 
linearly from about 1.3 GeV to approximately 0.1 GeV for proton energies between 6000 
and 7000 GeV. 

Another important experimental factor is the detector alignment. It is required that 
the detector stations will be able to measure the scattered proton trajectory elevation 
angles with precision of about 1 /irad. This implies a 10 /im precise alignment. It turned 
out that the 10 /im misalignment of the stations results an offset of the reconstructed 
proton energy. This offset has the largest value of about 5 GeV for protons of 6000 GeV 
energy and decreases to zero with proton energy increasing to 7000 GeV. 

5. Central Exclusive Production (CEP) 

The AFP detectors can be used to measure the exclusive central production of scalar 
J PC = ++ particles (for example the Higgs boson or some supersymmetrical particles). 
The central production can be viewed as a two stage process. In a first step each of 
the incident protons emits a color singlet object. Subsequently, these objects interact 
with each other giving a centrally produced system. The incoming protons remain intact, 
traverse the magnetic lattice of the machine inside the beam pipe and can be detected in 
the AFP detectors. The centrally produced system decays into the ATLAS main detector. 
Hence, this gives a unique possibility to measure all the particles belonging to the final 
state (a completely exclusive event measurement). 

Such events were simulated in a simplified way. In the generation the four momentum 
transfer, t, and the reduced proton energy loss, £, were distributed according to e~ bt with 
b = 6 GeV~ 2 and respectively. Later, the proton transport to the AFP220 detectors 
was simulated using the FPTrack calculations. 

In Figure M the geometrical acceptance for different masses of the centrally produced 
system for various distances between the detector edge and the beam centre is shown. As 
expected, the geometrical acceptance strongly depends on this distance and for a realistic 
distance of 3 mm (marked with the dashed line in Figure [6]) it varies between and 30% 
for masses changing from 300 to 800 GeV. 

Next, the mass of the centrally produced system was estimated using the detector mea- 
surements and the proton energy reconstruction described in section HI For the Central 
Exclusive Production process the produced system mass determination from the reduced 
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200 400 600 800 

Missing mass [GeV/c 2 ] 

Figure 6. The geometrical acceptance of detectors at 216 m away from the ATLAS 
Interaction Point for Central Exclusive Production process as a function of the produced 
particle mass. The solid lines depicts the acceptance for active detector region at the 2 
mm distance from the beam centre. The dashed, dotted and dash-dotted lines mark the 
acceptance curves for the 3, 4 and 5 mm distance, respectively. 

energy losses of both protons, £i and £ 2 , is possible via [[TU]: 

M x = yjs ■ 6 ■ 6, 

where s is the centre of mass energy squared. 

The mass reconstruction resolution as a function of the centrally produced system be- 
tween 300 and 800 GeV is shown in Figure [7] for the 3 mm distance between the detector 
and the beam centre. The impact of several experimental factors is also depicted in 
this figure. The mass reconstruction resolution, after an initial jump at the acceptance 
edge, very slowly increases from 5 to 8 GeV with increasing value of the produced sys- 
tem mass. The influence of the multiple scattering, the beam energy variation and the 
proton direction angular spread is small and below 2 GeV. In fact, the resolution value 
is dominated by two factors. First one is the detector spatial resolution which gives the 
contribution ranging between 2 and 6 GeV and which dominates for masses above 500 
GeV. The second one is the uncertainty on the xq coordinate of the interaction vertex, 
whose influence practically does not depend on the produced system mass and which is 
the most important factor for masses below 500 GeV. The field imperfections, estimated 
as described previously, have a small influence on the reconstructed mass resolution which 
is about 0.7 GeV at 300 GeV and saturates at the value of approximately 1 GeV at 500 
GeV. Also, in this case the field variation resulted the mass offset which is about 1 GeV 



Transport Simulation and Diffractive Event Reconstruction at the LHC 



11 



> 

a 

.2 
'-+-> 

o 

x 

CD 

f-H 



CD 
X 

a3 



all effects 
detector resolution 
xO unknown 



B field variation 
multiple scattering 
beam E and 6 smearing 



10 



300 400 500 600 700 800 

Missing mass [GeV/c 21 



Figure 7. The centrally produced particle mass resolution determined with outgoing 
protons as a function of the particle mass. The overall mass reconstruction resolution is 
marked with the thick solid line, the influence of: the detector spatial resolution - the 
thick dashed line, the vertex position in the transverse plane - the thick dash-dotted line, 
the multiple Coulomb scattering - the dotted line, the magnetic field variation - the solid 
line, the beam energy and the proton direction angular spreads - the dashed line. 



in the considered mass range. The impacts of the interaction vertex position and that 
of the detector misalignment, not shown in Figure are small and below 0.5%o of the 
produced mass value. The detector misalignment introduces the reconstructed mass shift 
which almost linearly increases with the produced mass value from about 2.5 GeV at 300 
GeV to 6.5 GeV at 800 GeV. 



6. Summary and Conclusions 

A parameterisation of the proton transport through the magnet lattice of the LHC was 
devised. This parameterisation has a simple functional form and enables fast and easy 
calculations. 

A proton energy unfolding procedure from the proton trajectory position measurements 
was prepared. This procedure allows the reconstruction of the scattered proton energy. 
The procedure was used to reconstruct the missing mass of the centrally produced scalar 
system. The missing mass reconstruction resolution weakly depends on the produced 
mass and reaches about 8 GeV at the mass value of 800 GeV. 

The proton energy unfolding procedure can be used for the first level triggering of the 
apparatus at the LHC environment. 
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